BCS superfluidity in ultracold gases with unequal atomic populations 
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We consider the existence of a BCS superfluid phase in ®Li due to the pairing of two hyperfine 
states with unequal number of atoms. We show that the domain of existence for this phase will 
be increased to a very large extent in the vicinity of the threshold for collapse. This is due to the 
presence of a new phase with anisotropic order parameter. This phase is induced by the anisotropic 
part of the scattering, linked to the indirect interaction between atoms through density fluctuation 
exchange. 

PACS numbers : 03.75.Fi, 32.80.Pj, 67.90.+Z, 74.20.Fg 



O 



o 



^ . 
i> ■ 
o ■ 
o : 
o . 



I 

O 

o 



X 
J3 



Much progress has been made recently toward the ex- 
perimental investigation of degenerate ultracold Fermi 
gases. Indeed the degenerate regime has already been 
reached in ^"K by sympathetic cooling of two hyperfine 
states 10] whereas the possibility of cooling a mixture of 
^Li and ^Li atoms in order to reach this range has been 
demonstrated |^ . In addition to a number of very inter- 
esting physical effects linked to phase space restriction 
10,^ one of the most fascinating prospect for experiment 
is the possibility [Q to reach the onset of Cooper pairs 
condensation and observe the resulting BCS superfluid, 
the equivalent of Bose Einstein condensates for Fermi 
gases. As we will see this phase is likely to present quite 
unusual features. 

In addition to the problems raised by the use of evapo- 
rative cooling for Fermi gases a possible stumbling block 
in the path to the BCS superfluid is the need to achieve 
a near equality |^,^ between the number of atoms in the 
two hyperfine states assumed to form pairs. Indeed in 
the various schemes under investigation there is no fast 
relaxation mechanism which equalizes these two popula- 
tions and the number of atoms in each state will be es- 
sentially conserved. Therefore in these ultracold atomic 
gases, one will necessarily have to deal with situations 
where the two populations of particles assumed to form 
pairs are not equal. Now the difference in chemical po- 
tential between these two populations acts as an effective 
field which tends to break pairs and the superfluid phase 
is destroyed when this difference is of order of the criti- 
cal temperature. Hence this pair breaking effect becomes 
more important a problem if the critical temperature is 
low and, since there is presently some uncertainty in 
the value of Tc, it is clearly of interest to investigate this 
question in detail. Moreover, independently of this possi- 
ble experimental problem, this difference in chemical po- 
tential is an additional control parameter in the system 
quite interesting to play with and it is worth studying its 
effect, all the more since it is not trivial. 

This problem of pairing with two unequal populations 
has already been considered a long time ago in the con- 
text of superconductivity. Here this is an applied mag- 



netic field which tends to make the two spin populations 
unequal. Usually the critical field is limited by orbital 
effects. However the question of the limitation of the su- 
perconducting phase, when these orbital effects are small, 
has been considered very early and it was pointed out by 
Clogston and Chandrasekhar that the standard 
BCS phase could, at most, resist to a difference in chem- 
ical potential between the two populations of the order 
of the critical temperature (the so-called paramagnetic 
limit). Not long after, Fulde and Ferrell P] (FF), and in- 
dependently Larkin and Ovchinnikov ]l(]| ] (LO), showed 
that pairing could somewhat adjust to the difference in 
chemical potential, instead of just resisting, by letting 
the pairs have a common momentum K instead of hav- 
ing it equal to zero as for a standard superconductor. 
However the effect was actually found rather small. In- 
deed at zero temperature the standard BCS phase goes 
to the normal state by a first order transition J7| when the 
chemical potential difference 2/2 is equal to v2Ao , where 
Aq = l.lQTc is the zero temperature gap for equal pop- 
ulations. The FFLO phase goes to the normal state by 
a second order transition for p, = 0.754Ao, which is not 
much beyond. Actually there is to date no undisputed 
observation of this phase in standard superconductors 
|]ll| , most likely because it is quite sensitive to impuri- 
ties. It is the purpose of the present paper to show that 
pairing can adjust even better than in the FFLO phase, 
and that this new phase can in some instances lead to 
quite a strong increase of the existence domain for the 
superfluid phase. 

One can view the appearance of the FFLO phase in the 
following way. When the chemical potentials for the two 
spin populations are different, it becomes more costly in 
terms of kinetic energy to form (k, — k) pairs in the stan- 
dard BCS way because one can not pick the two particles 
very near the Fermi surfaces, since these surfaces do not 
match due to their size difference. In the FFLO phase 
this problem is remedied by taking a nonzero total mo- 
mentum which amounts to shifting, in momentum space, 
one of the Fermi surface with respect to the other so that 
they almost match on some region. But naturally match- 
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ing gets worst on the opposite side. Nevertheless the total 
balance is barely favourable. Now one can think to im- 
prove this situation by making pairing stronger on the 
side where there is energy gain and weaker where there 
is energy loss. This means one looks for anisotropic pair- 
ing to find a lower energy ground state. Naturally this 
can not work if scattering is isotropic, as it will be es- 
sentially the case for a weakly interacting ultracold gas 
since p-wave scattering is negligible, and this will be in- 
deed the situation in the dilute regime where the cou- 
pling constant A = 2kp\a\/iT is small. However, even if 
the bare scattering is isotropic, the renormalized interac- 
tion is not because of the existence of the Fermi surface. 
This effect is at the basis of the Kohn and Luttinger 
^ paper, where they showed that, even with a repul- 
sive interaction. Cooper pairs would necessarily form in 
high angular momentum. For our purpose the case of 
^Li is of particular interest since the high density regime 
is bounded by an instability Q| occuring for A > 1 . In 
the vicinity of this instability the effective interaction will 
be strongly anisotropic, leading to a marked increase in 
the existence domain of the superfluid phase, compared 
to the FFLO phase. Note that most experiments are 
likely to be done in this range since it corresponds to 
higher critical temperature. We remark also that a sim- 
ilar phase with anisotropic order parameter should exist 
in clean superconductors with Pauli limited upper critical 
field and anisotropic interaction. 

Before going into the effect of anisotropic pairing, it is 
of interest to show that this is the best possible choice 
for pairing with unequal particle number. Indeed assume 
that we look for the most general kind of pairing, where 
k "f is paired with /(k) |. On one hand we want pairing 
to be stable against scattering, which means that scatter- 
ing must send this pair into another pair k' "f , /(k') |. 
On the other hand scattering conserves momentum which 
leads, for any k and k', to k -f- /(k) = k' -I- /(k') = K 
where K is a constant. This shows that /(k) = K — k, 
that is pairs with a total momentum K is the most gen- 
eral solution. Note that this argument assumes transla- 
tional invariance. This will not hold for a finite sample 
such as those obtained with trapped gases, which will in- 
duce surface effects. However these surface effects should 
be small if we do not want to have Tc drastically reduced 
(more precisely we want the pair size to be small 
compared to the sample size). So it is reasonable to ne- 
glect size effects in a first step. Now the common mo- 
mentum K produces a breaking of rotational symmetry, 
and the standard symmetry analysis leading to pairs with 
a given angular momentum (s-wave pairs in the case of 
^Li) is no longer valid. The general order parameter At 
will depend on the wavevector k. However we still have 
rotational invariance around K and we can classify the 
solutions by their angular momentum m. We will as- 
sume that the most stable pairing corresponds to m = 0, 
as it is likely to be so for a standard interaction. How- 



ever other values of m, corresponding to a breaking of 
the rotational symmetry around K, do not seem to be 
excluded from first principles and would be undoubtedly 
a very interesting situation. 

We will explore now quantitatively the possibility of- 
fered by an anisotropic order parameter. However our 
purpose is more to demonstrate the importance of the ef- 
fect than to perform an exact calculation. This last goal 
would require a perfect knowledge of the effective interac- 
tion, which is not available. We will consider specifically 
the case of ^Li and we restrict ourselves to the determina- 
tion of the T = critical difference in chemical potential 
2/2 = /i I — [ above which superfluidity disappears. 
Since we want to go in the high density regime, we need 
to have an expression for the effective interaction in this 
range. For this purpose we take the paramagnon model 
which we have already used for an evaluation of in the 
high density regime Actually we will not retain the 
full interaction of this model because most of the terms 
lead to a moderate anisotropy and they would produce 
a more complex calculation without adding much effect. 
The essential contribution of those terms is already in- 
cluded in the value of the critical temperature for 
equal populations. We will only retain the explicit at- 
tractive part coming from density fluctuations which pro- 
duces near the instability a strong contribution for low 
momentum transfer and in this way lead to a strongly 
anisotropic interaction. Moreover we will for simplicity 
perform a weak coupling calculation, omitting all the fre- 
quency dependence and taking the zero frequency value 
of the interaction. To be coherent with this weak cou- 
pling approach we will also omit self-energy effects. This 
leads to the total effective attractive interaction ^(k, k') 
H given by : 



^/nk,k') = A + i/^°(^ 



Axo(q) 



(1) 



with q = k — k'. Here Nj = mkp/2-K'^ is the density 
of states at the Fermi surface for equal population and 
Xo(q) is the reduced elementary bubble at zero frequency: 



2X0(q)^l + i(l-^)ln|±^ 

y 4 2-y 



(2) 



with y = q/kp. The first term in Eq.(l) is the direct 
term and the last one is the indirect interaction due to 
density fluctuation exchange. 

The difference in chemical potential 2/2 has the effect 
of shifting the kinetic energies of the particle (measured 
from chemical potential) by ±/2. For the pair propaga- 
tor this is just equivalent to shift the frequency uo 
by /2. Similarly taking the momenta of the pair mem- 
bers to be ±k -I- K/2 instead of ±k produces a shift by 
±k.K/2m, leading to an overall shift in frequency by 
fik = P- — k.K/2m. At finite temperature this means re- 
placing the Matsubara frequency ujn = 7rT(2n -I- 1) by 
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LUn — ijlk- On the other hand we expect the modification 
of the efi'ective interaction caused by ji to be small, of 
order fi/Ep, and we can neglect it. This leads us finally 
to the following gap equation at finite temperature : 



dk 



^2 



(1 + 



1 Axo(q) 



AI + (LJn-ifik)^ ' 2 l-Axo(q) 



:) (3) 



where is the kinetic energy measured from the Fermi 
surface for p, = 0. Since we are looking for the criti- 
cal jl and expect a second ordre phase transition as for 
the FFLO phase, we let Ak go to zero in the denom- 
inator. Then, in order to get rid of the cut-off in the 
frequency summation, it is convenient to make use of 
7rrX;(l/|wn|) = + ln(rO/T), where is the critical 
temperature for /I = 0, K = and in the absence of the 
indirect term in the interaction. When we specialize to 
T — > 0, the frequency summation can be performed by 
^T^sgnK)[(l/K-i/2k)-l/t^n] =ln(0.88T/|/ik|). In 
this way we obtain the following equation for the critical 
diff'erence in chemical potential jl : 



A, 



^A,(l + Aln^)(l + 1^^^ 
47r ''^ ^ 2\flkr 2 l-Axo(q) 



) (4) 



where Ag is the zero temperature gap for /I = 0, K = 
and in the absence of the indirect term in the interac- 
tion. When the indirect interaction is not present, Ak is 
independent of k. Then the angular integration is eas- 
ily performed, leading to ln(Ao/2/l) = (l/2)/o(x), with 
X = Kkp / i'^mp,). Here fo{x) = j^^du ln|l + xu\ which 
is negative and minimum for x w 1.200 (x = cotha;) 
leading to the standard FFLO result Ji — 0.754Ao . 

We consider now the effect of the indirect interaction. 
Since both k and k' are on the Fermi surface, q goes from 
to 2k F and xo(q) goes from 1 to 1/2 . Although it is 
not a problem to use numerically the exact expression 
for xo(q), one can sec that xo(q) = 1 — a + a cos 6, with 
a = 0.25 and q = 2kF iim(0/2) is a quite good approx- 
imation. Indeed it gives properly the two limits q = 
and q = 2kF- Moreover it corresponds to keep only the 
first two terms in a Legendre polynomials expansion with 
slightly modified coefficients (the exact I = coefficient 
is (1 + 21n2)/3 = 0.755 and the exact coefficient of cos^^ 
is 0.232 , the higher order terms being fairly small). This 
approximation allows to perform the azimuthal integra- 
tion analytically, and since our calculation is anyway a 
model calculation it is quite reasonable to make this sim- 
plification, despite its slight inaccuracy in the vicinity of 
61 = 0. 

However, before proceeding to the solution of the 
resulting equation, it is interesting to consider first a 
slightly simplified version which allows to carry out 
the calculation completely explicitely. We merely re- 
place the interaction term by the first two terms of 



its Legendre polynomial expansion Vq + Vi cos 9, with 
Vo = 0.5 + 0.25(aA)-ilog[l + 2aX/{l - A)] and Vi = 
-1.5(aA)-i +0.75(1 - A + aA)(aA)-2 log[l + 2aA/(l - A)]. 
In this case the solution of Eq.(4) has the form Ak = 
1 + Si cos a with cos a = k.K. Then L = ln(2/i/Ao) is 
solution of the second order equation (L— 1/A+l/ XVq + 
0.5/o)(L -1/X + S/XVi + 1.5/2) = 0.75/?. We have set 
fn{x) = f^^du In |1 + a;u|. fi{x) and f2{x) are both 
found to be largest for x in the range 1.1 - 1.2, and the 
largest L is also found for the same range of x for all 
values of A. 
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FIG. 1. Pull line : location of the second order transition 
/Lt/Ao to the normal state from Eq.(4) as a function of A. 
Dashed line : approximate solution with Legendre polyno- 
mial expansion. Long-dashed line : location of the first order 
transition p,/Ao = 1/V2. 



Naturally the indirect interaction produces a trivial 
effect, namely the renormalization of the coupling con- 
stant A into AVo, producing a corresponding change of 
the critical temperature and of the gap for equal popula- 
tion, which goes from Aq to Aq = Aq exp(l/A — I/AVq). 
This effect is simply found by looking at the solution Lq 
for X — 0, that is without FFLO phase, since in this 
case only the isotropic part of the interaction is relevant. 
In the present case this is given by = 1/A — 1/AVb. 
We arc only interested in i — Lq which gives the in- 
crease in the superfiuid domain due to the existence of 
our FFLO phase, in units of Aq. In Fig.l we have plot- 
ted as the dashed fine /x/Aq. We see that, roughly for 
A < 0.6, there is essentially no change with respect to 
the standard FFLO result. And indeed the gap remains 
essentially isotropic in this range. On the other hand, for 
A > 0.6, p, increases rapidly and for A = 0.9 it is almost 
twice the standard FFLO result. It is also quite interest- 
ing to consider the anisotropy linked to ^i. As soon as 
p starts to grow with respect to the standard FFLO re- 
sult, the order parameter gets a sizeable anisotropy. For 
A = 0.94 we obtain a node on the Fermi surface, and for 
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larger A a change of sign over the Fermi surface with a 
nodal line. 




FIG. 2. Angular dependence of Ak for A = 0.6 , 0.7 , 0.8, 
0.9 , 0.95 and 0.99 as a function of the angle a between k and 
K (A increases with A for a = 0). Inset : anisotropy of the 
order parameter as a function of A. 

Let us turn now to the results of the numerical solu- 
tion of Eq.(4) . They are given as the full line on Fig.l 
for /i/Ao as a function of A (as above we have taken into 
account the renormalization of Aq into Aq by calculating 
L — Lo). For small and intermediate A the result is es- 
sentially identical to the result of our simplified version. 
However when A approaches 1, p, increases more rapidly. 
This is easy to understand by looking at Eq.(4). In this 
regime the effective interaction becomes very large for 
small q and the log term favors wavevectors nearly along 
K so it is better to have the gap function Ak peaked 
for wavevector along K (see Fig. 2). This can be seen 
more explicitely by looking at the simplified form taken 
by this equation when A — > 1. Taking into account the 
numerical evidence that in this limit the optimal p, is 
obtained for x = I, and making the change of variable 
k.K ~ u(l — X)/a, we obtain for F{u) = Ak: 



F{v) 



1 

'4^ 



^ i7t ^ M + \ogu 

^l + 2{u + v) + {u~vY 



(5) 



with M — L + log((l — A)/a) — 1. This equation can 
be solved numerically. However the solution we are look- 
ing for behaves approximately as exp(— /iu). Inserting 
this form into Eq.(5) and requiring that Fiv) is zero 
for large v gives log/i = M — C (for a — 0.25), where 
C is the Euler constant. On the other hand requiring 
F{0) = 1 and making an asymptotic evaluation of the 
resulting integral gives /i = exp(— \/2) ~ 0.24 which 
is (surprisingly) not much different from the numerical 
result n ~ 0.3. This leads finally to the asymptotic 
evaluation L w — 1 — log(l — A), which, together with 
Lo ~ 1 — l/log(0.5/(l — A)), is in reasonable agreement 
with our direct solution of Eq.(4) found in Fig.l . Natu- 



rally the divergence of ^ itself is not to be taken seriously 
since our calculation requires jl/Ep to be small anyway. 

Coming back to Fig.l , we see that beyond A « 0.8 
the results from the full Eq.(4) get larger than those of 
the simplified equation. Hence we find a very large in- 
crease of the domain for our FFLO phase. When we 
take into account |6) that the critical temperature itself 
will increase rapidly in this range due to the indirect in- 
teraction, we see that looking in this region seems quite 
promising experimentally since the overall domain for the 
superfluid phase will be much increased. On the other 
hand the physical properties of this phase will be to a 
large extent quite different from those found for equal 
population. A first reason is that this phase is gapless, 
just as the standard FFLO phase |^,|l^. Next we find an 
important order parameter anisotropy. Indeed we have 
plotted in Fig. 2, for various values of A, the angular de- 
pendence of Ak with respect to K. As we have men- 
tionned already, Ak gets more concentrated along the K 
direction when A increases. The value of its minimum 
compared to its maximum is also plotted in Fig. 2 . We 
see that even for moderate values of A such as A « 0.4 the 
anisotropy is quite sizeable. However when A is further 
increased the anisotropy becomes ultimately huge to the 
point that Ak ~ for a very large fraction of the Fermi 
surface. This is a highly unconventional situation and 
certainly quite unique among BCS superfluids. Another 
more well-known feature will further complicate the mat- 
ter. We have in our system a degeneracy with respect to 
the direction of K, which will be in general lifted by a 
texture leading to a spatial inhomogeneity, that is yet 
another symmetry breaking, as investigated by LO [Tc| ] 
and more recently in Ref. ||ll|. Together with experi- 
mental inhomogeneity due to the trap this will lead to a 
remarkably complex physical situation. 

In conclusion we have shown that in ®Li the indirect 
interaction due to density fluctuations exchange will lead 
to the appearance of a new BCS phase with anisotropic 
order parameter. Near the instability threshold, this re- 
sults in a large increase of the superfluid domain as a 
function of the difference between the atom numbers in 
the two hyperfine states forming the Cooper pairs. 
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